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The preceding Comment HI] makes a few claims on the in- 
serted figure of Fig. 5 in our Letter 101 (denoted as FIG- 
URE hereafter). The FIGURE displays the variations of the 
superdiffusion exponent a versus the screening parameter k 
for different effective coupling strength Toff ISfl . The authors 
in HI] first claim that our observations in this figure "resulted 
from an incorrect account of the coupling strength" and in- 
stead propose a new definition of so-called relative coupling 
strength F™' = T/Tc, where Fc is F at melting point. How- 
ever, it is not hard to prove that F'''^' HI] is essentially equiv- 
alent to our Feff li: F^' = T/T, = [F/(k)]/[F,/(k)] = 
rofr(K)/Foff,c, where /(k) is a scaling formula invented by 
Kalman et al. Since Fcff,c is a constant in Jst], F''''' and 
Foff are related with a constant coefficient Fcff,c- Second, the 
authors claim that the formula that we adopted in our Letter 
to calculate Foff is only valid when k < 3. Unfortunately, 
there were no explicit or implicit statements in Ref. that 
this formula is only valid for k < 3. This might be a newly- 
discovered issue for this formula, as we understand. In that 
case, it only affects the last two data points at k = 3.5 and 4.0 
in the FIGURE of whereas the claim of "incorrect account 
of the coupling strength" seems exaggerated. 

In the Comment the authors also re-produce the FIG- 
URE ]|2] in terms of F"^"' and make comparison between our 
results and theirs. We suppose that in their calculation ex- 
actly the same parameters as ours ^ are used, except for the 
fitting range, which is Updi £ [100 320] in iflf]. In our Let- 
ter ||2[], we use different fitting ranges depending on the sys- 
tem states (Details of fitting technique will be presented else- 
where.). In particular, for the data shown in the FIGURE, we 
use ujpdt G [220 320] and ujpdt e [100 320]; the former gives 
the values of a and the difference of the two fits gives the un- 
certainties of a. The reason for doing so is clearly shown in 
Fig [U in which one observes that the tail of MSD /< is not a 
straight line (in log-log scale). Instead its slope changes with 
time, so different fitting ranges will certainly result in differ- 
ent values of a (This tendency is more significant for smaller 
F and was also observed by Donko et al. 10] recently.). To 
mimimize this effect, we choose a fitting range as close to the 
end of the tail as possible to approximate the long time asymp- 
totic behavior of MSD / 1 (Note that a is defined by the slope 
of the MSD/t in long time limit.), but a wider range to esti- 
mate its uncertainty. This can explain some of the differences 
between our results and theirs for k < 3.0. The dramatic 
differences for k = 3.5 and 4.0 are indeed due to the differ- 



ent system states. We have performed new calculations with 
given F'°' in Comment and the results are shown in Table 
I. However, we still observe a detectable difference, especially 
forF'''' = 0.75: our results are substantially larger than theirs. 
The task of resolving all these differences is much beyond the 
scope of this Reply. Nevertheless, the detailed algorithms for 
our simulation was available in Ref. ][5j and the source code 
for our Letter ^ will be provided upon request. 

In addition, the authors of the Comment ]|H claim that the 
dependence of a on k is "regular and systematic" based on 
their observations. It is a very interesting extension to the re- 
sults of ]|2] and could shed new light on the effect of interac- 
tion stiffness on superdiffusion. However, the given explana- 
tion is not satisfactory and there is no any relation between 
their observations and the explanation. 
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FIG. 1: MSD{t)/t (normalized by cjpda^) for k = 4.0, and F = 
338 (corresponding to F'^' = 0.075 of H). Dash-lines are fits of 
asymptotic behaviors: the upper one is fitted in the range ujpdt G 
[220 320] but is plotted to the full range to ease the comparison, and 
the lower one in the range LOp^t G [100 320]. a is the slope of the 
fits. 
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TABLE I: a for k — 3.5 and 4.0. Numerators and denominators 
are values of a fitted in the ranges ujpdt € [100 320] and ujpdt £ 
[220 320], respectively 



«:\r''°' 


0.075 


0.375 


0.75 


3.5 


0.18/0.154 


0.166/0.149 


0.069/0.062 


4.0 


0.203/0.184 


0.164/0.156 


0.058/0.054 
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